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We consider a localized electronic spin controlled by a circularly polarized optical beam and an external 
magnetic field. When the frequency of the beam is tuned near an optical resonance with a continuum of higher 
energy states, effective magnetic fields are induced on the two-level system via the Inverse Faraday Effect. 

We explore the process in which the frequency of the beam is made linearly time-dependent so that it sweeps 
through the optical resonance, starting and ending at the values far away from it. In addition to changes of spin 
states, Kramers-Kronig relations guarantee that a localized electron can also escape into a continuum of states. 

We argue that probabilities of transitions between different possible electronic states after such a sweep of the 
optical frequency can be found exactly regardless the shape of the resonance. We also discuss extension of our 
results to multistate systems. 
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I. INTRODUCTION 

Light interaction with a magnetic system is sensitive to optical polarization. For example, absorption of a circularly polarized 
light in a semiconductor depends on the magnetic moment of conduction electrons ID. In addition to such dissipative effects, 
a polarized optical beam can experience more subtle magneto-optical effects, such as rotation of the linear polarization axis of 
a beam, called the Faraday Rotation. It is frequently used in optical spectroscopy la. A related effect is the shift of the phase 
of the circularly polarized light passing through a medium with spin polarized electrons 0. Microscopically, such effects are 
described by virtual dissipationless excitations of the magnetic degrees of freedom. Despite seemingly different microscopic 
origin of dissipative and nondissipative magneto-optical effects, causality of the magnetic medium response to the interactions 
with light imposes constraints between these effect. Such constraints are known as the Kramer-Kronig relations lH. 

As magnetic media influence the polarized light, there is a feedback from this interaction on the dynamics of the electron 
spin polarization. Moreover, both dissipative changes of light intensity and dissipationless phase shifts of the beam induce 
counterpart effects on electronic spins, such as dissipative optically induced electronic spin polarization in semiconductors 11 
and dissipationless appearance of an effective, optically induced, magnetic fields acting on electronic spins. The latter effect is 
generally known as the Inverse Faraday Effect 0. It has attracted considerable interest recently due to applications to purely 
optical ultrafast control of magnetization 0. The held of research called the circuit QED Q aims in engineering such light 
interactions with discrete quantum states in order to realize desired entanglement and dynamics of solid state qubits. 

Coupling of quantum spin systems to polarized light can be strongly enhanced when the frequency of light is tuned near an 
optical resonance with excited states of electrons. However, unlike the application of time-dependent fields, Kramers-Kronig 
relations guarantee that in addition to effective time-dependent fields, optical pulses inevitably induce dissipative transitions JEl 
that break coherence or even remove electrons, e.g., from a confining potential of quantum dots. Hence, in order to optimize the 
optical qubit control, it is important to develop the theory of qubit interactions with time-depenent optical pulses at near resonant 
conditions. 

In this article, we derive an exact nonperturbative result, which is akin to the famous Landau-Zener formula 0. The latter 
is the solution of the model that describes interaction of a two-level system with a linearly time-dependent magnetic field. 
This model is used frequently in applications of nonstationary quantum mechanics, e.g., for characterization of parameters of 
nanoscale systems inoHni and developing quantum control strategies MM- 

In our case, instead of a linearly time-dependent magnetic field, we consider a situation when a spin system interacts with a 
polarized optical beam such that the frequency of the beam changes linearly with time and sweeps through an optical resonance. 
As in the case of the Landau-Zener model, when the frequency of the beam is far away from the resonance, interactions of the 
beam with the two-level system are negligible. Only during a finite time interval, when the frequency of the beam is close to the 
resonant frequency, can the beam induce additional dynamics of a two-level system. Hence, one can define a scattering problem; 
given the state of the spin system and t — oo, what is the state of the spin at f -t-oo? Using the approach pioneered by 

Landau M and some of the recent developments M, we demonstrate that one can derive simple exact analytical expressions 
for transition probabilities between spin states in this scattering process. 


II. THE MODEL 


We consider a situation shown in Eig. [T] Circularly polarized light interacts with spin up (but not down) state of a two-level 
system. The frequency of the light w is tuned to be near the resonant frequency ujq of a transition to a continuum of high energy 
states. If such a transition happens irreversibly, electron is lost in the continuum. In addition to such irreversible transitions, 
electron can produce virtual coherent transitions that merely renormalize the energy of the spin up state. 

One can describe the interaction of a spin-1/2 with the light plus a constant external magnetic field h by the Schrodinger 
equation 


with an effective Hamiltonian; 




H = 


x{^) 0 
0 0 


+ h • O’, 


( 1 ) 


( 2 ) 


where x(a;) is the response function that describes renormalization of the energy of the spin up sate. This function is complex 
valued due to the possibility of irreversible excitations. Let us explicitly introduce the real part of this function Xi(‘^) its 
imaginary part X 2 (w); 


X(w) = xi(w) + iX2(w)- 


( 3 ) 
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FIG. 1. (Color online) An electron with spin 1/2 is initially in a superposition of two spin states with energies split by a magnetic field. A 
circularly polarized beam with frequency lj interacts with this electron. The optical frequency uj is tuned near the resonant frequency tJo that 
corresponds to the excitation of an electron into a continuum of much higher energy states. We assume that if an irreversible transition in the 
continuum happens, the electron has negligible probability to return to the initial doublet during the time of experiment. Frequency of the 
beam is assumed to be linearly time-dependent cu — tuo = pt, sweeping from large negative to large positive t with a constant rate /3 > 0. 
Circularly polarized field has nonzero matrix elements between the states of the continuum and the spin “up” state. Spin “down” state does 
not interact with this field. 


Karamers-Kronig constraints relate them as 


Xi(w) = 

-V . 


(4) 
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X2(w) = - 


r 
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where the symbol V near the integral means that only the principle part of this integral is considered. 

In what follows, we assume an almost arbitrary form of the function x(w) except that its dissipative part is localized in the 
frequency domain, so that the integral 

k = J duj'x2iuj') ( 6 ) 


is finite. The constant k characterizes the strength of the resonance. 

In appendix A we work out a specific case of a broad resonance with a weak beam intensity, for which expressions for Xi{^) 
and X 2 (w) can be obtained explicitly, and Kramers-Rronig relations can be verified. It is important to note here that the Kramers- 
Kronig relations can be justified physically in the context of Eq. Q only if dissipative electronic transitions are truly irreversible, 
which means that being excited to the continuum, the electron must leave the system. Such an experimental situation can be 
found in the literature, e.g., on spin decoherence in InGaAs quantum dots ini. On the other hand, if the excited electron can 
return to the quantum dot during the time of experiment, one should use a more complicated density matrix approach, to which 
our discussion does not apply. In what follows, we will focus on mathematical properties of the dynamics described by equations 

Q-®- 

The relation between causality, analyticity and Karamers-Kronig relations, which is discussed thoroughly, e.g., in Q, pre¬ 
dicts that any function x(a;) that satisfies Eqs. ffl-® must be analytic function of w in the upper-half of the complex plane 
(Re(a;),Im(a;)). Setting w —/ oo, we can disregard u}^n denominator of Q and find that 

k 

Xi(w) ~ -, |w-wo|»7, (7) 

w — Wo 

where 7 is a characteristic width of the resonance in the frequency space, and fc is a constant given by Eq. •S'- Thus we find that 
the real part of x(w) decays as ~ l/w at w —)■ cxd. 

Eor example, let the absorption be described by a Lorentzian: 


X 2 (w) = - 


kx 


(w-wo)2+72 


( 8 ) 


where wq 7 is the optical frequency of the resonance. The minus sign in (j^ is needed to guarantee that the Schrodingier 
equation Q leads to a decay of the amplitude. Substitituting into Q we find: 


Xi(w) = 


k{ui — ujq) 

(w - wo)2 +7^ 


(9) 
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FIG. 2. (Color Online) The real and imaginary parts of x(^) for ^ Lorentzian (Eqs. ^ - l |10^ ) and a Guassian Xg(w) (Eqs. 0 

and ( |12t ) cases with the same strengths k — 1 and widths 7 = 1 of the resonances. The Gaussian absorption decays much faster than the 
Lorentzian one. 


Equations ([^ and (|^ can be written more compactly as 


x(w) 


k 

(w — Wo) + ^7 


( 10 ) 


which makes it clear that the pole of ( [TOl i is at the point w = wq — 17 , i.e. it is placed at the lower complex half-plane. 
Generally, we assume that x(w) can have multiple poles or other singularities in the lower half-plane. For example, consider a 
Gaussian-broadened shape of absorption: 


then 


X2(w) 



Xi(w) 



^ erfi 


LJ — CJq 

V2x 


( 11 ) 


( 12 ) 


where erfi(a;) = dy is the imaginary error function. For any large w in the upper half of the complex plane, including 

the real axis, we have: Xi{^) ~ k/{uj — uJo) at|w —wqI ^ 7 , i.e. it satisfies our physical assumptions. Such a function x(w) has 
no singularities at finite values of w but it has exponentially diverging singularity at w —>^ 00 in the lower complex half-plane. In 
Fig. I^we compare real and imaginary parts of the Lorentzian and Gaussian response functions. In particular, this figure shows 
that while imaginary parts have strongly different behavior at large w, the real parts of x(^) h^ve similar asymptotics. 

The form of x(^) can possibly be very complex. For example, due to the linearity of Kramers-Kronig relations, any function 
that is a sum of several above mentioned functions x(w) taken with different values of constant parameters also satisfies our 
constraints. 


III. SWEEP OE THE BEAM EREQUENCY 

Imagine now the situation when the frequency of the laser beam is made linearly time-dependent: 

LU — LUq = pt. (13) 

Alt ^ ± 00 , the imaginary part of x{^) strongly suppressed, and the real part behaves as ~ 1/t, which is still negligible in 
comparison to the constant coupling to the external magnetic field. Therefore, if a system is prepared in one of the eigenstates 11) 
or 12) of the Hamiltonian (j^ with xi^^) 0’ transitions between such states would not happen at f —± 00 , and the scattering 

problem is well defined for a sweep of frequencies throughout the resonance. 

Suppose the system is initially (at t —> — 00 ) in the eigenstate 11) of the Hamiltonian, and we can measure the state of electron 
in the quantum dot after the sweep, at f > -Foo. There are 3 possibilities: First, electron can simply escape from the dot into the 
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FIG. 3. (Color online) The time evolution contour along the real axis (black) can be continuously deformed into a semicircle with radius 
i? —>■ oo laying in the upper complex half-plane (green) without affecting initial conditions at real values of f —>■ — oo. Since the function x{t) 
is analytic in the upper-half plane, this deformation does not encounter the singularities of the Hamiltoninan so that the scattering matrix does 
not change during the time contour deformations. Analogous continuation to the lower half-plane would be impossible because of singularities 

ofx(f). 


continuum. Second, it can stay in the dot and remain in the same state |1). Finally, it can stay but flip the spin into the state |2). 
All those processes have some probabilities to be realized. Let Pi^j be the probability that the system, starting in the state |f), 
will end up in the state \ j), where i,j = 1, 2. Probabilities of all possible outcomes are described by the following 2x2 matrix: 

P = ( \ (14) 

Pi^2 P2^2 )■ ^ ’ 

Note that this matrix is sufficient to determine the probability to escape, e.g. 

Pl^il = 1 ~ Pl^l — Pl^2, ^2->.0 = 1 ~ ^2-).l — P2^2, (15) 

where 0 denotes the state with electron in the continuum. 

Despite the simplicity of the sweep almost arbitrary complexity of the function x(a;(f)) makes the problem to determine 
the scattering probability matrix ( [T4| ) look extremely complex. For example, we know that one cannot solve a problem of a 
2-state system driven by arbitrary time-dependent fields analytically and exactly. Nevertheless, in this article we show that 
no matter how complex the function x(w) can be, as far as it satisfies the above mentioned constraints, we can provide exact 
formulas for three out of four elements of the matrix ( [l4| ). 

IV. SOLUTION OF THE MODEL 

We now solve the model using a contour deformation method described in M- This method takes roots in the Landau’s 
derivation of the Landau-Zener (LZ) formula ca. First, since we are looking at transitions between eigenstates of the magnetic 
held, it is better to switch to the basis of these states. Let 6 be the angle between the magnetic held and the beam axis. Then the 
Flamiltonian reads: 


H = (i + cos(6»)(T2 + sm{e)a^) , (16) 

where e is the parameter that describes the energy splitting by the external magnetic held and 1 is the unit 2 x 2 matrix; uj{t) is 
given by Without a loss of generality, we will set e > 0. 

Since the Hamiltonian ( [T6| ) has no singularities in the upper plane, the solution of the Schrodinger equation is analytic in this 
region. Hence we can deform the time contour continuously so that it has the form of a semicircle C with radius i? 3> wq, 7 , e, 
as shown in Fig.[^ This contour connects the same initial and final points at real values of t ^ ± 00 , so the solution along this 
contour can be used to solve the desired scattering problem. 

Along the contour C, the function x(a;) behaves as 


X(w)|c = ^+o{l/R), 


( 17 ) 
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where k is the real constant defined previously. 

Let, initially, the system be in the state |1) with energy e/2. Following ifTfil . we can say that since the off-diagonal elements 
of the Hamiltonian along the contour C become negligibly small, we can apply adiabatic approximation everywhere along 
the contour C, which corresponds in our case to keeping only contributions to the adiabatic energy that decay not faster than 
^ 1/R at i? —^ (X). In particular, this means that we can disregard the off-diagonal terms of the Hamiltonian. The adiabatic 
approximation leads then to the solution for the amplitude to remain in the state 11) along the contour C; 


■01 (r) = exp 



A:(l-I-cos(0))' 
^ I3t 



(18) 


An analogous solution is obtained for the amplitude to remain in the state |2) if the spin system is prepared in this state at 

t —)■ —oo: 


-02 (t) = exp -i 


oo+zO ^ 


fc(l — cos( 6 *)) 




dt 


(19) 


We note that although the form of x(w(f)) can be very complicated, Eq. (18 1 is valid as long as Eq. (|T^ is satisfied, since the 
arguments leading to that equation only depend on asymptotic behavior of the Hamiltonian at i? —oo. 

As for the amplitude of the transition from one spin state to the other state, trimming off-diagonal elements of the Hamiltonian 
for the evolution along the contour C leads to the prediction of a zero transition probability between states |1) and |2). However, 
according to Ref ca, this expectation is valid only for the transition from the state | 2 ), which has the lower adiabatic energy 
along the contour C, to the state |1). In terminology of ifTa . the transition |2) —^ |1) is “counterintuitive”. The discussion 
of exact suppression of counterintuitive transitions and the reasons why the adiabatic approximation fails in our case for the 
transition from | 1 ) into | 2 ) would be completely analogous to the discussion of the no-go theorem in ifTfilfTSl . 

By performing integration in ([T8]l-([T9| along the full contour C, and taking into account the no-go theorem, we obtain the 
following scattering matrix elements: 

S'i 2 = 0, 5-22 = (20) 


Taking the absolute value squared of the scattering matrix elements, we find the form of the transition probability matrix: 

g-7rfc(l-|-cos(6())//3 Q 

g-irfe(l-cos(6())//3 

where X indicates that this element of the transition probability matrix cannot be determined by the above described procedure. 

Equation ( |2T] l is the main result of this article. It shows that three out of four independent transition probabilities are described 
by simple formulas. In Eig.Qa-d) we show numerical results for time dependence of transition probabilities in the models with 
Gaussian and “Double-Lorentzian”. In particular, they demonstrate that at intermediate time moments the transition from the 
state 12 ) into the state | 1 ) can be nonzero but it approaches zero value asymptotically at f -|-oo, i.e. the zero value of the 
corresponding matrix element in ( [2T| i is not trivial. In Eig. j^e) and (f) we compare predictions of Eq. ( [ST) with results of our 
direct numerical simulations. 

It is likely that a simple exact analytical expression for the remaining unknown element of the transition probability matrix 
0 does not exist. In order to shed light on its behavior, we identified one fully solvable model with the Hamiltonian of 
the type ( [T 6 | ), which complete solution we worked out in Appendix B. The model corresponds to the Lorentzian shape of the 
dissipative part of the function x(^) (Eqs. ([^ - ([TOll). The corresponding Schrodinger equation can be mapped to the confluent 
hypergeometric equation, which solution has well-understood asymptotics. The result for the non-universal component then 
reads: 

X = — e-^fc(i-cose)/^j^ (-22) 

Equation explicitly shows that X depends on both e and 7 . Eigure shows that the matrix element denoted by X is not 
universal. It depends sensitively on all parameters of the model, in particular on the choice of the functional form of x(a;). 
Nevertheless, other scattering probabilities are independent of e (as far as it is positive), as well as numerous details that can 
describe the function x(^)- 




V. EXTENSION TO MULTISTATE SYSTEMS 

Our results are straightforward to generalize to multistate systems linearly driven through a resonance lfT9ll20l . Eor example, 
consider a three-level system (possibly a spin-1 multiplet of an electron in a quantum dot molecule) with a Hamiltonian 

/ x(w) 0 0 \ 

0 0 0 

\ 0 0 0 J 


H = 


(23) 
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(a) Guassian 



(c) Double-Lorentzian 



(e) Gaussian 



(b) Guassian 



t 

(d) Doubie-Lorentzian 



t 

(f) Doubie-Lorentzian 



FIG. 4. (Color Online) Numerical tests of analytical results. The dots correspond to direct numerical solutions of the Schrddinger equation 
0. and the solid curves correspond to theoretical predictions in Eq. \2l\ . In all cases, we make the choice of parameters /3 = e = 1. (a) Time 
dependence of the four transition probabilities for the Guassian case. Choice of parameters: A: = 1, 7 = 0.1, 0 = 57 r/ 12 . The time-interval 
for simulations is from t = —50 to 50 and the time step is 0.0002. Only the curves in the interval from t = —10 to 10 are shown. Each 
probability undergoes a sharp change at f = 0 and then oscillates before saturation, (b) Time dependence of the probability P 2 ^i zoomed 
in near f = 0. Note that although it approaches zero asymptotically at f —>■ + 00 , it is nonzero at intermediate time moments, (c) Same as 
(a) but for the “Doubie-Lorentzian” case: x(^) = \i ui-u}i+i-y ui-u} 2 +i-y ^' where we choose cui = —u >2 = 1. (d) Same as (b) but for the 
Doubie-Lorentzian case. Note the existence of two peaks, (e) Transition probabilities as functions of the strength k for the Gaussian case. 
Choice of parameters: 7 = 1, 0 = tt/S. The time-interval for simulations is from t = —500 to 500 and the time step is 0.002. The numerics 
show perfect agreement with theoretical predictions, (f) same as (e) for the Doubie-Lorentzian case. 
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FIG. 5. (Color Online) Numerical results (the dots) for the probabilities Pi _>2 (or X) for three different cases with the same strength k = 1: 
Lorentzian, Gaussian and “Double-Lorentzian”. The choice of other parameters k P = 'y = t = 1, 9 = j, and uji = —a ;2 = 1 for 
double-Lorentzian. The time-interval for simulations is from t = —500 to 500 and the time step is 0.002. For the Lorentzian, solid curve 
represents the theoretical prediction of Eq. \22\. 


For simplicity, we assumed here that the light interacts with only one of the states. In the basis of eigenstates of the Hamiltonian 
with x(a;) = 0 the Hamiltonian matrix reads; 


H = 


e 0 0 \ 
0 0 0 
0 0 -e ) 


x(^) 

4 


1 -V2 1 \ 

-V2 2 -V2 

1 -V2 1 / 


(24) 


Repeating the same arguments as in the two-level case, and using the no-go theorem from m, we find the form of the transition 
probability matrix: 


/ 0 0 \ 

P = Xi^2 0 . (25) 

V Xi^3 X2^3 J 

There are three unknown elements in ( [25] l but all other transition probabilities are determined exactly regardless of the details of 
x(a;). The numerical check for this model, in support of Eq. ( [25] ), is presented in Fig.j^a). Figurej^b) shows numerical results 
of the three unknown transition probabilities. Equation ( |25| ) can be readily generalized to an arbitrary n-level system interacting 
via a resonant transition. 

We would like to note that there are potential applications of this multistate generalization to the field of ultracold atomic 
gases interacting via the Feshbach resonance. The latter corresponds to a coherent creation of molecules from pairs of single 
atoms. The strength of atomic interactions via such a resonance can be tuned by a magnetic field and has a typical ~ 1 /{B — Bq) 
dependence on the external field. There already have been numerous experiments performed to study the effect of a linear sweep 
of the magnetic field through the Feshbach resonance in order to create a molecular condensate ED. Certainly, direct application 
of our results to such systems will need additional justification because creation of a molecule may not be possible to describe 
as a dissipative transition. However, phenomenologically, if we identify the discrete states with quantum states of an atomic 
gas, our results predict, in particular, that the sweep through the Feshbach resonance cannot leave the atomic part of the gas in 
an excited state. Moreover, we can obtain the formula for the probability to remain in the initial state after the sweep through 
the resonance. If the atomic gas is initially in the ground state, then combining such predictions we can derive the probability 
of not creating any molecule after the sweep through the Feshbach resonance e.t.c.. In fact, the possibility of such results is in 
agreement with some of the existing theoretical publications on this topic casi. 


VI. CONCLUSION 


We considered a localized electronic spin controlled by a circularly polarized beam and an external magnetic field. We 
assumed that the frequency of the beam is made linearly time-dependent so that it sweeps through the optical resonance, starting 
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FIG. 6. (Color Online) Numerical results for a three-level system with Hamiltonian Eq. where x is taken to be Lorentzian. The 

dots correspond to direct numerical solutions and the solid curves correspond to theoretical predictions in Eq. Choice of parameters: 
/3 = 7 = e = 1. The time-interval for simulations is from t = —500 to 500 and the time step is 0.002. (a) Analytically determined transition 
probabilities as functions of the strength k. Here Pi^i and coincide. Pa-ji and P 3_»2 are not shown but both of them are tested to 

be zero, (b) The three non-universal transition probabilities as functions of k. Pi _»2 and P 2-^3 coincide. The inset shows Pi-s-a which is 
significantly smaller than the other two probabilities. 


and ending at the values far away from it. Causality of the spin response to the optical interactions imposes the Kramers- 
Kronig relations on the response function and make it possible for the electron to escape into a continuum of states. Using a 
contour deformation method, we found exact probabilities of transitions between different possible electronic states after such a 
sweep of the optical frequency. These probabilities do not depend on the shape of the resonance but rather they depend on its 
integrated power k of the resonance and the relative angle between the optical beam and the applied external magnetic field. This 
simplicity can be used in practice, e.g. to estimate the parameter k experimentally. We showed that this result can be generalized 
to multistate systems which can be useful, e.g., for studies of transitions through the Feshbach resonance in ultra cold atomic 
gases. We believe that our theoretical predictions can be verified experimentally in spin systems controlled by optical beams, 
e.g., in NV centers |0 or in electron doped InGaAs self-assembled quantum dots iflTll . 

Acknowledgment. The work was carried out under the auspices of the National Nuclear Security Administration of the U.S. 
Department of Energy at Los Alamos National Laboratory under Contract No. DE-AC52-06NA25396. Authors also thank the 
support from the LORD program at LANL. 


APPENDIX A: LOCALIZED ELECTRON COUPLED TO A CONTINUUM OE CONDUCTION STATES 


As a physical example, consider an electron state |l/2) localized in a quantum dot of a semiconductor interacting with a 
circularly polarized beam. This light couples the state |l/2) to the trion exciton with the total spin 3/2. We assume that this 
excited state is not stable so that if it is created, one electron quickly escapes into the conduction band leaving zero charge in 
the quantum dot, which makes this quantum dot optically inactive during the time of experiment. Therefore, it is convenient to 
write the exciton state in the basis of a continuum of states of conduction electrons, |i), where index i runs throughout all Bloch 
states. The Hamiltonian of interaction of the localized electron and the continuum can be written then in the form: 

4ot=u;i/2|l/2)(l/2|+^L;,|/)(z|+ |^e*‘"‘^p,|z)(l/2| + /i.c.^ , (26) 

where are the matrix elements of the dipole operator between localized electron and conduction band states. Such couplings 
are also proportional to the beam field amplitude. Here we set the energy of the state | — 1/2) to zero. Since it does not interact 
with the given circalar polarization, this state drops out from ( |26] l. As shown in Lig. we assume that the gap between states 
with energies Ei and uJi/2 is in the optical spectrum so that uj ^ Ei — uJij2 ^ \Ej — Ei\. We will also assume that the couplings 
gi are weak in comparison to the broadening of the exciton state. 

Our goal is to derive the effective Hamiltonian describing the localized spin. Lor this, first, we redefine the phase of the band 
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states: |i) e This makes the Hamiltonian time independent: 

Htot = W 1 / 2 II/ 2 ) (1/2| + ^^{Ei — a;)|i)(t| + y^gi\i){\/2\ + h.c. 

i \ i 

We use then the stationary perturbation theory to derive the renormalization of the energy of the spin state |l/2): 



(27) 


LUi/2 —>■ LOi/2 — 'P 


dEp{E) 


IgjEr 

E — UJ — U!ij2^ 


(28) 


where we switched from the sum over discrete index i to the integral over energies E of delocalized states, and introduced the 
density of states p{E) and the energy-dependent coupling strength g{E). 

Equation (j2^ includes only virtual transitions to the exited states. However, there are inevitable dissipative transitions, which 
rate is given by the golden rule: 


C^IV’1/2P 

dt 


= -27r j dEpiEME)\^d{E-Lo-Lo,f2), 

and which can be introduced in the Schrodinger equation by adding a complex energy to the state |l/2): 


uJi /2 uJi /2 - iTT j dEp[E)\g{E)\^5 [E - uj - ^ 1 / 2 ) . 


Combining Eqs. (28 1 and (30 1 we hnd uJi /2 Wi /2 + x{^)’ where 

\9m^ 


X(uj) = -'P f dEp{E) -jTT f dEp{E)\g{E)\'^S {E - u; - UJ 1 / 2 ) , 

J £/ — W — U!i/2 J 


(29) 


(30) 


(31) 


where the first and the second terms in ( p^ can be identihed with, respectively, the dissipationless Xi{^) dissipative X 2 (w) 
parts of the response function x(a;). To show that they satisfy the Kramers-Kronig relations, it is sufficient to consider a special 
case p{E) = 6{E — Eq) (because Kramers-Kronig relations are linear, and an arbitrary density p{E) can be represented as a 
continuum of such delta-functions). In this case 


x(w) = |g(^^o)|' 


1 


- inS (Eq -w- a;i/ 2 ) 


Eq-U - LOi/2 

It is then straightforward to verify that relations (j^ hold for corresponding Xi(^) X 2 (w). One can also rewrite 

uj — h/Q -\- oj-]_j2 H“ 

so that x(<^) clearly has only a single pole in the lower half of the complex plane. 


(32) 


APPENDIX B: TRANSITION PROBABILITIES FOR THE LORENTZIAN SHAPE OF THE RESONANCE 


Consider the Hamiltonian •mi) with a Lorentzian shape of X 2 (w). Under the linear sweep of the frequency. 


X{t) = k/{l3t + ix). 


We will assume that e > 0. Let us write the state vector as ip = {(pi,ip 2 )’^- The change of variables t ^ t — ij//3 eliminates 
the parameter 7 from the evolution equation (but moves it to the boundary conditions). Eliminating ipi from the Schrodinger 
equation (0, we get a 2nd order differential equation for ip 2 - 


tip 2 + 




ek cos 9 

w 


+ 



(P2 = 0 . 


(33) 


This equation can be transformed to the confluent hypergeometric equation: 


zc"{z) -I- (r — z)c'(z) — sc{z) = 0, 


(34) 
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where z = —iet, c{t) = e r = 1 + and s = ^ which has a solution 

ip 2 {t) = e~^ [AM (s, r, r) + B17(s, r, r)], (35) 

where t = —iet + M{s,r,T) and U{s,r,T) are the Kummer’s and Tricomi’s functions, respectively, and A and B are 
constants to be fixed by initial conditions. The solution for the other state amplitude reads: 

0 T T 

(pi(t) = — tan -e~"^ {A[M{s, r, r) + -M(s + 1, r + 1, r)] + B[U{s, r, r) — tU{s + 1, r + 1, r)]}. (36) 

The asymptotic behaviors of the functions M(s, r, t) and U (s, r, t) at large |t| read: 

M{s,r,T) ^ 

1 (s) i (r — s) 

U{s,r,T) (37) 

where the cut of arg r is taken to be the negative real axis and arg(—1) = tt. They are used to determine, under specified initial 
conditions at f —)• —oo, the constants A and B. Under the initial conditions, (^i(—oo) = 0 and |(^ 2 (—oo)| = 1, we obtain A = 0 
and |i 7 | = e'^'''/( 2 / 3 )—irfe(i-cose)/ 4 / 3 ^ asymptotic behaviors at f —>• oo read: 1 ^ 1 ( 00 ) = 0 and |(/ 32 (c!o)| = 

Squares of moduli of (pi(oo) and (p 2 {oo) then give the transition probabilities from the second state: 

P2^1 = 0 , 

P 2^2 = (38) 


Under the initial conditions, |(pi(—oo)| = 1 and (p 2 (— 00 ) = 0, we obtain |A| = |cos||e ^ 
B = — A g—7rfc(i-cos e)/2/3^ rpjjg ti-ansition probabilities from the first state then read: 


p _ —7rfc(l+cos6>)//3 

^ 1—>-1 — C , 


(39) 


The results ( [38l l and ( |39| l agree with Eq. (|2T| obtained via the contour deformation method. In addition we have obtained the 
non-universal element X = Pi^ 2 - 
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